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ABSTRACT 

The projected structures and integrated properties of the Andromeda I, II, III, V, 
VI, VII and Cetus dwarf spheroidal galaxies are analysed based upon resolved counts 
of red giant branch stars. The observations were taken as part of the Isaac Newton 
Telescope Wide Field Survey of M31 and its environs. For each object, we have de- 
rived isopleth maps, surface brightness profiles, intensity-weighted centres, position 
angles, ellipticities, tidal radii, core radii, concentration parameters, exponential scale 
lengths, Plummer scale lengths, half-light radii, absolute magnitudes and central sur- 
face brightnesses. Our analysis probes to larger radius and fainter surface brightnesses 
than most previous studies and as a result we find that the galaxies are generally 
larger and brighter than has previously been recognised. In particular, the luminosity 
of Andromeda V is found to be consistent with the higher metallicity value which 
has been derived for it. We find that exponential and Plummer profiles provide ade- 
quate fits to the surface brightness profiles, although the more general King models 
provide the best formal fits. Andromeda I shows strong evidence of tidal disruption 
and S-shaped tidal tails are clearly visible. On the other hand, Cetus does not show 
any evidence of tidal truncation, let alone disruption, which is perhaps unsurprising 
given its isolated location. Andromeda II shows compelling evidence of a large excess 
of stars at small radius and suggests that this galaxy consists of a secondary core com- 
ponent, in analogy with recent results for Sculptor and Sextans. Comparing the M31 
dwarf spheroidal population with the Galactic population, we find that the scale radii 
of the M31 population are larger than those for the Galactic population by at least 
a factor of two, for all absolute magnitudes. This difference is either due to environ- 
mental factors or orbital properties, suggesting that the ensemble average tidal field 
experienced by the M31 dwarf sphcroidals is weaker than experienced by the Galactic 
dwarf sphcroidals. Wc find that the two populations are offset from one another in 
the central surface brightness - luminosity relation, which is probably related to this 
difference in their scale sizes. Finally, we find that the M31 dwarf spheroidals show 
the same correlation with distance-from-host as shown by the Galactic population, 
such that dwarf spheroidals with a higher central surface brightness are found further 
from their host. This again suggests that environment plays a significant role in dwarf 
galaxy evolution, and requires detailed modelling to explain the origin of this result. 

Key words: Local Group - galaxies: general - galaxies: dwarf - galaxies: fundamental 
parameters - galaxies: structure - galaxies: interactions 



1 INTRODUCTION 

Most of our detailed knowledge on the structure of dwarf 
spheroidal (dSph) galaxies comes from observations of the 
nine 1 dSp hs which make up p a rt of t he Galactic satellite sys- 
tem. Ilrwin fc Hatzidimitrioul il995T) (hereafter IH95) anal- 

1 IWillman et alj 120051) have recently discovered a new Galactic 
companion in Ursa Major, which is not included here 



ysed the structure of the Galactic dSphs, except Sagittar- 
ius, by mapping their resolved star counts from photographic 
plates. They found that the stellar profiles of the dSphs were 
generally well described by a single component King or ex- 
ponential model. A generic feature was an excess of stars at 
large radii, over and above that expected from the best-fit 
King tidal model, which have generally been interpreted as 
evidence for tidal disturbance. A King profile is not a phys- 
ically motivated model for dSphs since their relaxation time 
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is of order a Hubble time and their stellar velocity distribu- 
tion may deviate significantly from Maxwellian. Any inter- 
pretation of the structure of a dSph based exclusively on a 
King model fit must therefore be treated with caution. Nev- 
ertheless, such analyses do yield a useful parameterisation, 
and in particular the equivalent half-light radius can read- 
ily be compared with alternative parameterisations derived 
from Plummer or exponential profiles. The King model tidal 
radii (r t ) of the Galactic dSphs are of order 1 kpc, and their 
average half-light radius (ri) is ~ 170 pc, although Fornax 
is over twice this size. More recent work on the radial pro- 
files of Galactic dSphs from wide field CCD cameras (Ursa 
Minor: [Maiewski et alJl200Ct iMartmez-Delgado et alJl200lt 
Draco: loderiHrcn^n^t^uTl200ll : Carina~ |Palm^^t^lTll2003i 
generally confirm the results from earlier photographic stud- 
ies. 

Some recent studies have revealed that the structures 
of dSp hs are more complex than had previously been as- 
sumed. [Sarbecj^^ii] (l200lf) looked for spatial variations in 
the colour of the RGB and the horizontal branch of several 
Local Group dSphs and found evidence for hor izontal branch 
popul ation gradients in many of their systems. iTolstov et alJ 
(2004) have shown that Sculptor has a spatially, chemically 
and kinematically distinct second componen t in the form 
of a ce ntrally concentrated core. Additionally. iKlevna et ail 
(2004) have shown that the Sextans dSph has a kinemati- 
cally cold core, and that Ursa Minor has a distinct concen- 
tration of stars offset fro m its geometric cen tre which has 
distinct kinematics (IH95: Fklevna et alJl20o3) . 

Given these recent findings in the Galactic dSph sys- 
tem, it is timely to look at the next closest dSph system 
to our own, that of M31. Approximately 16 satellites make 
up this subsystem, of which seven - Andromeda I, II, III, 
V, VI, VII, and IX - are classified as dSphs. Andromeda I, 
II and III were discovered by Sidney van den Bergh in his 
pioneer ing survey for Local G r oup g alaxies in the early 
1970's llvan den Berghl Il972blal Il974h . together with An- 
dromeda IV, a n object which was la ter shown to be a back- 
ground galaxy llFerguson et alJl2000j) . Andromeda V and VI 
were discovered some years later bv lArmandroff et alJ ( 1998, 
ll99ST) from a detailed analysis of the digitised Second Palo- 
mar Sky Survey. At the sa me time, Andromeda VI was dis- 
covere d independently by iKarachentsev fc Karachentseval 
( 1999) along with with Andromeda VII. An additional ob- 
ject, Andromeda VIII llMorrison et al.l|2003f) . has recently 
been postulated to be associated with an over-density of 
planetary nebulae i n the giant stella r stream visible in the 
south-east of M31 dlbata et alJl200ll) . although the nature 
of this over-density remains unclear. An dromeda IX was re- 
cently discovered bv lZucker et alJ i2004f) and is the faintest 
satellite of M31 yet known, with M v ~ —8.3. 

In comparison to the dSphs which orbit the Galaxy, 
relatively little is known in detail about the stellar popu- 
lations of the Andromeda dSphs. From the colour of their 
red giant branches, it would ap pear that all of them ar e 
metal poor [Fe/H] < -1.5 dex f eg. IMcConnachie et ai]l2005l) . 
Deep HST fields for Andro meda I, H an d III reaching be - 
low the horizontal branch <Da Costa et alll99d.l200CU2002l) 
show extended epochs of star formation and variations in 
horizontal branch morphology, suggesting the star forma- 
tion histories have been notably different. In addition, An- 
dromeda I is observed to display a gradient in its horizontal 



branch morphology, such that their are more blue horizon- 
tal branch stars located at larger radius from the centre of 
the dwarf llDa Costa et alJll996ft . Some evidence for AGB 
compone nts have also been seen i n th ese systems (most re - 
cently bv lKerschbaum et aljl2004l and lHarbeck et al-lEooll . 
although a strong intermediate age component similar to 
some of the Galactic dSphs is generally lacking. 

As a compliment to the deeper pointed HST observa- 
tions, we have obtained wide field Johnston V (V) and 
Gunn i (i r ) photometry for the majority of the members 
of the M31 subgroup using the Wide Field Camera (WFC) 
on the 2.5m Isaac Newton Telescope (INT). This is a four- 
chip EEV 4 K x 2K CCD mosai c camera with a ~ 0.29D 
field of view dWalton et alJl200ll) . With typical exposures of 
~ 1000 s in each passband, this photometry is deep enough 
to observe the top few magnitudes of the red giant branch 
in each system and has already been used to derive a homo- 
geneous set of distance and metallicity estima tes for each of 
these galaxies llMcConnachie et all20oll200Sll Here, we use 
the same data to analyse the structural properties of each of 
the Andromeda dSphs, as well as the isolated dSph in Ce- 
tus as a comparison. In common with Tucana, this galaxy 
is one of only two dSphs not found as a satellite to a large 
galaxy in the Local Gro up. Our results for A ndromeda IX 
are presented elsewhere dChapman et al.ll2005l) . Our overall 
technique and methodology is similar to that adopted by 
IH95 for the Galactic dSphs, ins ofar as we base most o f our 
analysis on reso lved star counts. TCaldwell et alJ dl992l) and 
ICaldwell lll999l) have performed a similar analysis for the 
Andromeda dSphs based upon the integrated light, and we 
will compare our results to these later (Section 4.1). 

In Section 2, we derive contour maps, radial profiles 
and associated structural parameters for Andromeda I, II, 
III, V, VI, VII and Cetus based upon resolved star counts. 
The integrated luminosities, central surface brightnesses and 
related quantities are derived in Section 3. We postpone dis- 
cussion of all the results until Section 4, which also compares 
the M31 dSph population to the Galactic dSph population. 
Section 5 summarises and concludes. 



2 STAR COUNTS AND STRUCTURE 
2.1 Preliminaries 

The reader is referred to [ Ferguson et alJ ([2002) and 
IMcConnachie et alJ (120041 . 12003) for information on our ob- 
serving strategy and data reduction process. We do, how- 
ever, recap here the morphological photometric classification 
procedure for each detected source, as this is central to our 
subsequent analysis. 

Object Classification: Objects are classified indepen- 
dently in each passband based on their overall morpholog- 
ical properties, specifically their ellipticity as derived from 
intensity-weighted second moments and the curve-of-growth 
of their flux distribution (llrwin et alJl2004T) . Measures from 
these are combined to produce a normalised N(0,1) statis- 
tic which compares the likeness of each object to the well- 
defined stellar locus visible on each frame. Stellar objects are 
generally chosen to lie within 2 or 3 a of this locus depending 
on the desired trade-off between completeness and contam- 
ination from non-stellar objects. Contamination takes the 
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form of a small number of spurious images which are es- 
sentially eliminated by requiring objects to be detected on 
both passbands. Additionally, at faint magnitudes (within 
~ 1 — 2 mags of the frame limit depending on seeing), dis- 
tant compact galaxies may also be incorrectly identified as 
stellar images. This is particularly true in bluer passbands. 
For this specific study, completeness is more important than 
mild non-stellar contamination. We therefore use all objects 
identified in the INT WFC that lie within 3 a of the stellar 
locus in the i'-band, and require only that the object has 
been detected in the V'-bwad within 1 arcsec of the i'-band 
position, to reduce contamination from background galaxies. 

Background corrections: A good estimate of the back- 
ground contamination and its uncertainty is vital to accu- 
rately determining the extent and profile of each dSph. The 
term "background" describes both Galactic foreground stars 
and distant compact galaxies. These are minimised in the 
first instance by using simple cuts in colour - magnitude 
space, designed to isolate the RGB and remove those stars 
whose colour and magnitude make membership of the dSph 
unlikely. This is illustrated for Andromeda II in Figure Q 
A horizontal cut in i' magnitude, coincident with the de- 
rived location of the TRGB, removes objects too bright to 
be RGB stars in the dSph. Two cuts on the blue and red 
side of the RGB also removes objects whose colour makes 
membership of the dSph unlikely. This process removes a 
significant fraction of the Galactic foreground. 

Although the influence of compact galaxies is signifi- 
cantly reduced by the morphological classification/selection 
scheme, it cannot be completely removed at faint magni- 
tudes. A statistical correction must therefore be applied to 
remove the remaining contamination. This is estimated by 
excising the dSph from a pixelated map of the stellar spatial 
distribution. An intensity distribution is constructed from 
the remaining pixels and a sigma-clipped least-squares fit 
of a Gaussian profile is performed on this "background" 
distribution to calculate the position of the peak (ie. the 
background level) and associated errors in the usual way. 
As we do not know the full extent of the dSph a priori, the 
correct subtraction of the dwarf component requires several 
attempts before we are left with the easily identifiable back- 
ground component to perform the fit to. We do not require 
to simultaneously fit the dwarf profile plus background be- 
cause the background is relatively flat over the area of the 
field, with no evidence of strong differential extinction for 
any of the dwarfs which would affect this measurement. For 
large dSphs which cover the majority of the field of view (An- 
dromeda II, VII and Cetus), the background estimate will 
suffer from contamination from the outer regions of these 
objects. This effect is generally small, as the number den- 
sity of stars belonging to the dSph at large radius is much 
smaller than the background contamination. However, this 
effect will lead to a slight overestimate of the background 
by a few percent, resulting in a marginal steepening of the 
outer parts of the derived profile. The level of this effect 
is small enough to be covered by the derived uncertainties 
on the scale-radii, and will always act to underestimate the 
actual values of these parameters. 

Crowding corrections: In the central regions of the dSphs, 
the higher density of images causes some individual sources 
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Figure 1. A colour magnitude diagram of Andromeda II in Lan- 
dolt V and I. The dashed lines represent the cuts that we use 
in order to reduce background contamination. The horizontal cut 
is at the location of the tip of the red giant branch derived in 
McConnachie et al. (2004, 2005). The other two cuts are placed 
so as to isolate the red giant branch loci as shown. 



to remain unresolved. This can be approxi mately corrected 
for by using the crowding correction of llrwin fc Trimble) 
(1984). 



/' - / 1 



(1) 



where / is the observed number density of all images (not 
just stellar), /' is the actual number density and A' is the 
typical area of the image, where the radius of the image 
is closely approximated by the seeing. For the INT WFC 
observations, the typical seeing is ~ I arcsec. The deriva- 
ti on and validity o f Equ ation [T] is given in the Appendix 
of llrwin fc Trimble! lll984t) and ignores second-order effects 
such as the actual shape of the luminosity function. For the 
dSphs, the typical background stellar densities are ~ 2 - 3 
stars arcmin~ 2 (after photometric cuts have been applied) 
and are unchanged by the crowding correction. The crowd- 
ing correction increases the stellar counts in the central re- 
gions of the dSphs by typically 5 — 10 %. 



2.2 Isopleths and dwarf geometry 

A field of 0°5 x 0°5 containing the dSph is used for each 
isopleth map. For Andromeda II, V, VI, VII and Cetus, 
this corresponds closely to the INT WFC pointing, and the 
north-west corner is missing in each case due to the geometry 
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of the CCDs. For Andromeda I and III, multiple pointings of 
the camera are combined to form a mosaic. For each object, 
the field is divided into ~ 100 x 100 pixels (~ 18arcsecs res- 
olution), and the crowding correction is applied based upon 
the source density in each cell. The isopleth maps of the stel- 
lar sources are constructed, smoothed and displayed in Fig- 
ure!^ For the two fields constructed from multiple INT WFC 
pointings, faint-end threshold magnitudes have been used to 
account for varying incompleteness levels between fields due 
to different observing conditions. Contour levels in Figure [5] 
are set such that each subsequent increment is a factor n 
greater than the previous increment, where the first incre- 
ment is 2 a above the background, n is different for each sys- 
tem and is given in the caption of Figure [5] This allows the 
same number of contours to be used for each dwarf galaxy, 
making visual comparison of their morphologies easier. The 
2 a contour also picks out the peaks of noise features as well 
as the dSph, and gives a general indication of the quality of 
the background subtraction. 

For each isophote, we calculate the centre of gravity 
(x , y ), position angle, 9 (measured east from north), and 
eccentricity, e, by using the intensity weighted moments. 
Thus, 
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where / (xt, yt) is the intensity of the i th pixel and a xx , o yy 
and <j xy are the intensity-weighted second moments within 
each isophotes. The eccentricity is related to the ellipticity, 
e = 1 — b/a, via 



6=1 



1 + e 



(3) 



The "average" parameters for each dSph are then estimated 
from the values derived for each isophote and the uncertainty 
estimated from the variation in the parameters as a function 
of isophotal threshold. This procedure gives estimates of 8 
and e that are independent of any parameterisation, and 
allows a more robust estimate of their uncertainties. Effects 
such as isophote rotation are likewise able to be quantified 
for those systems where such an effect is taking place. The 
results of this analysis are tabulated in Table 



2.3 Radial profiles and associated parameters 

The stellar radial profile for each galaxy is constructed by 
measuring the stellar number density in elliptical annuli with 
the average ellipticity and position angle derived above and 
applying the background correction. Assuming fixed param- 
eters for each galaxy gives a robust estimate of the profile 
suitable for comparison with other galaxies and models. The 
width of the elliptical annuli are set by requiring a minimum 
signal-to-noise threshold in each bin. The results are plot- 
ted in logarithmic form in Figure^] The error bars take into 
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Figure 3. - continued. 



account the Poisson error in the counts and the uncertainty 
in the background estimate. 

The expected form of the stellar density distribution in 
dSphs is unknown, although several profiles have been sug- 
gested. It is not our intention to compare the radial distribu- 
tions in Figure |3] to every model which has been proposed, 
and we have instead adopted the simplest, most intuitive 
and most commonly used. In particular, King profiles are 
used and the simple, empirical form of these is given by 



fK = A 



(4) 



\ [1 + (r/r c ) 2 ] 5 [1 + (r t /r c ) 2 ] 5 j 

jKindll962fi . A is a scaling parameter, r the radius from 
the centre of the system, r c the core radius and rt the tidal 
radius. For the range of core concentrations considered here 
(c = log 10 (r t /r c ) < 1.5), these pro files closely m atch the 
physically-motivated form derived in iKind l)l96ft) for glob- 
ular clusters, wh ich is based on a lowered Gaussian distri- 
bution function jBinnev fc Tremaindll987l) . King profiles - 
both the empirical and physically-motivated forms - provide 
a tractable family of models with intuitive par ameters that 
have been fitted extensively to dwarf galaxies llHodgelll966l : 
lEskridgdll988allbT: IH95). We therefore adopt them as one 
of the models to be compared to the radial distributions in 
Figure [3] although we note that there is not necessarily a 
physically-motivated reason to expect King profiles to be a 
realistic form of the radial profile for dwarf galaxies, and 
interpretations based on this assumption should be treated 

wit h caution. 

iFaber fc Linl il983l) advocate exponential profiles to de- 
scribe the projected surface density distribution of dSphs, 
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Shape Parameters 





HA (x ) 


Dec (y a ) 


"DA O 1 O \ 
PA = U { ) 


e = 1 — b/a 


Andromeda I 


00n45m40.3s 


+38 02 21 


22 ± 15 


0.22 ± 0.04 


Andromeda II 


01hl6m28.3s 


+33°25'42" 


34 ± 6 


0.20 ± 0.08 


Andromeda III 


00h35m31.1s 


+36°30'07" 


136 ± 3 


0.52 ± 0.02 


Andromeda V 


01hl0ml7.0s 


+47°37'46" 


32 ± 10 


0.18 ± 0.05 


Andromeda VI 


23h51m46.9s 


+24°34'57" 


163 ± 3 


0.41 ± 0.03 


Andromeda VII 


23h26m33.5s 


+50° 40' 48" 


94 ± 8 


0.13 ± 0.04 


Cetus 


00h26ml0.5s 


-11°02'32" 


63 ± 3 


0.33 ± 0.06 



Table 1. Intensity- weighted centres and geometric parameters for the dSphs. Errors are estimated from the variation in the parameters 
as a function of isophotal threshold. Centres are estimated to be accurate to within ±7 arcsecs in both directions. 



fs = B exp 



(•») 



where B is a scaling parameter and r e is the effective radius. 
These models require one parameter less than King models 
and qualitativel y often provide a s good a fit as the more 
complex family. iRead fc Gilmorel <l2004t) have shown that 
exponential profiles are a generic phase of dSph evolution if 
they undergo rapid mass loss at early times. r e is a useful 
parameter and we additionally compare this class of model 
to the stellar radial distributions. 

Finally, Plummer models are frequently used in N-body 



iFont e t al. 2004; La w et al.ll2005ft . They are of the form 

f P = c b2 



(b 2 + r 2 ) z 



(6) 



where C is a scaling parameter and b is the Plummer core 
radius. Due to their common use in simulations, we adopt 
this as the final model to compare with the dSph profiles. 

These three models are fitted to each of the stellar ra- 
dial profiles using a least-squares minimisation technique. 
The best-fitting King, exponential and Plummer models are 
overlaid on the logarithmic profiles in Figure[3]as the solid, 
dashed and dot-dashed curves respectively. The derived pa- 
rameters corresponding to each of these models is listed 
for each galaxy in Tables [5] and [3] along with the associ- 
ated uncertainties and the formal value of the reduced x 2 
statistic. The concentration parameter for the King profile is 
also given. The scale-radii correspond to the geometric mean 
along the two axes of the dwarf. The surface brightness scale 
on the right vertical axis of Figure [3] has been calculated by 
normalising the radial profile to the integrated flux measure- 
ments described in the following section. 



3 INTEGRATED PHOTOMETRY 

We have directly estimated the central surface brightness 
and integrated luminosity of the dSphs from the integrated 
flux distribution of each galaxy. With suitable image pro- 
cessing, the effects of Galactic foreground stars and random 



noise can be reduced to manageable levels, enabling these 
integrated parameters to be measured in a simple manner. 

The processing procedure is relatively straightforward. 
First, the existing object catalogues are used to define a 
bright foreground star component at least 0.5 - 1 mags above 
the TRGB, to allow for the potential presence of AGB stars. 
A circular aperture is then excised around each foreground 
star and the flux within this aperture is set to the local 
sky level, interpolated from a whole-frame background map. 
The size of this aperture is the maximum of four times the 
catalogue-recorded area of the bright star at the detection 
isophote, or a diameter four times the derived FWHM see- 
ing. Each frame is then re-binned on a 3 x 3 grid to ef- 
fectively create 1 arcsec pixels. The binned image is then 
further smoothed using a 2D Gaussian filter with a FWHM 
of 5 arcsecs. 

The result of this procedure is to produce a coarsely 
sampled smooth image containing both the resolved and 
unresolved light contribution from the dSph. The central 
surface brightness can then be trivially measured by deriv- 
ing the radial profile, here defined as the median flux value 
within elliptical annuli. Finally, large elliptical apertures are 
placed over the dSph and several comparison regions to es- 
timate the background-corrected integrated flux from the 
dwarf and the reference regions. The variation in the flux 
from the multiple comparison measures gives a good indica- 
tion of the flux error, which is, of course, dominated by sys- 
tematic fluctuations rather than by random noise. To mit- 
igate the effect of random residual foreground stellar halos 
and scattered light from bright stars just outside the field of 
view, the elliptical apertures are chosen to correspond to the 
derived value of ri for the dSph. The estimated total flux is 

2 

then scaled to allow for this correction. For Andromeda VII, 
a quarter-light radius is used instead as there is a significant 
Galactic nebulosity in this field which makes estimates based 
on larger apertures unreliable. As with the resolved compo- 
nent, some of the comparison regions may contain light from 
the dSph. However, this effect is negligible in comparison to 
the systematic fluctuation. 

The results for each galaxy are listed in Table |I| The 
integrated apparent magnitude has not been corrected for 



6 McConnachie & Irwin 




10 -10 10 -10 

S (arcmins) S (arcmins) 

Figure 2. Isopleth maps for the six dSph satellites of M31 (panels (a) - (f)), and the isolated dSph in Cetus (g), in standard coordinate 
(£, rj) projections. Each map shows a similar 0°5 X 0°5 area of sky. A large range in structural properties is clearly evident for these galaxies. 
The first contour of each map is 2 a above background. Contour increments then increase by a factor n for each subsequent contour: (a) 
Andromeda I: n = 1.15 (b) Andromeda II: n = 1.3 (c) Andromeda III: n = 1.15 (d) Andromeda V: n = 1.15 (e) Andromeda VI n = 1.25 
(f) Andromeda VII: n = 1.33 (g) Cetus: n = 1.25. 



extinction, but all other qua ntities are extin c tion-c orrected 
using the values gi ven by ISchlegel et al J (Il998t) . tabu- 
lated in Table 1 of iMcConnachie et alJ J2005I) . The dis- 
tance measurements and a ssociated uncertainties derived in 
IMcConnachie et alJ i2005l) have also been used. V is the 
directly measured central surface brightness for the dSphs. 
So is the derived value of the central surface brightness ob- 
tained by normalising the King profile such that the inte- 
gral under the profile is equal to the integrated magnitude 
for the dSph. For Andromeda II, a two-component profile 
has been used instead of the King profile (see Section [4. 3H . 
The uncertainty on the measurement of S is dominated by 
the uncertainty on the absolute magnitude for the dSph, 
and the uncertainty on the directly measured central sur- 
face brightness V is estimated to be ~ 0.1 — 0.2 mags. There 
is reasonable agreement between the values for V and S , 



which highlights the accuracy and consistency of the inte- 
grated magnitude measurement, the radial profiles, and the 
surface brightnesses. 

TableHJalso lists the derived values of ri of each dSph, 

2 

obtained by integration of the appropriate King profile. 
This is with the exception of Andromeda II, where a two- 
component profile was used, and Cetus, where ri has been 

2 

derived by integration of the exponential profile. The lack 
of an obvious truncation radius in the Cetus data results in 
a value for ri derived from the King profile, which is twice 
as large as from either the Plummer or exponential profiles, 
and which is probably unreliable. The half-light radius of a 
Plummer profile is equal to b, and the half-light radius of an 
exponential profile is 1.68 r e . Generally, these are consistent 
with the tabulated value for ri measured from the King 
model fit. 
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S (arcmins) 

Figure 2. - continued. 





x 2 


rc (') 


n (') 


King Profile 
r c (kpc) 


rt (kpc) 


c = lo Sio (rt/r c ) 


Andromeda I 


1.01 


2.7 ± 0.3 


10.4 ± 0.9 


0.58 ± 0.06 


2.3 ± 0.2 


0.59 ± 0.06 


Andromeda II 


1.84 


5.2 ± 0.2 


22.0 ± 1.0 


0.99 ± 0.04 


4.2 ± 0.2 


0.63 ± 0.03 


Andromeda III 


0.89 


1.3 ± 0.2 


7.2 ± 1.2 


0.29 ± 0.04 


1.5 ± 0.3 


0.74 ± 0.10 


Andromeda V 


1.02 


1.2 ± 0.2 


5.3 ± 1.0 


0.28 ± 0.04 


1.2 ± 0.2 


0.63 ± 0.11 


Andromeda VI 


0.96 


2.1 ± 0.2 


6.2 ± 0.4 


0.48 ± 0.06 


1.4 ± 0.1 


0.46 ± 0.06 


Andromeda VII 


0.91 


2.0 ± 0.1 


19.3 ± 1.6 


0.45 ± 0.02 


4.3 ± 0.4 


0.98 ± 0.04 


Cetus 


1.00 


1.3 ± 0.1 


32.0 ± 6.5 


0.29 ± 0.02 


7.1 ± 1.5 


1.40 ± 0.10 



Table 2. Details of the best-fitting King (Equation^ profiles shown in Figurc|3]as the solid curves. Each scale radius is the geometric 
mean for the dwarf. The distance moduli and associated uncertainties derived in McConnachie et al. (2005) have been used to transform 
units. 
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1 2 

radius (arcmin) 




2 5 
radius (arcmin) 



Figure 3. Log-log plots of the stellar radial profiles of the dSph galaxies. Error bars take into account the Poisson uncertainty in the 
counts and the uncertainty in the background estimate. The solid curves represent the best King profile fit to the data (Equation l4l 
Table |2"1 . the dashed lines correspond to the best exponential fit (Equation [F] Table 131 and the dot-dashed lines represent the best 
Plummer model fit (Equation |§] Table l3l . The surface brightness scale on the right vertical axis has been calculated by normalising the 
data to the integrated flux measurement in Section 131 



Exponential Profile 



Plummer Profile 





x 2 


r f 


(') 


r e 


(kpc) 


x 2 


6(') 


b ( kpc) 


Andromeda I 


1.36 


1.72 


± 0.06 


0.38 


± 0.02 


1.22 


3.12 ± 0.12 


0.67 ± 0.03 


Andromeda II 


1.66 


3.53 


± 0.06 


0.67 


± 0.01 


2.14 


6.44 ± 0.10 


1.23 ± 0.02 


Andromeda III 


1.00 


1.00 


± 0.05 


0.24 


± 0.01 


0.76 


1.82 ± 0.10 


0.40 ± 0.03 


Andromeda V 


1.25 


0.86 


± 0.05 


0.20 


± 0.01 


0.92 


1.56 ± 0.08 


0.35 ± 0.02 


Andromeda VI 


2.20 


1.20 


± 0.04 


0.27 


± 0.01 


1.56 


2.15 ± 0.08 


0.49 ± 0.02 


Andromeda VII 


1.14 


2.00 


± 0.04 


0.44 


± 0.01 


1.18 


3.47 ± 0.07 


0.77 ± 0.02 


Cctus 


2.06 


1.59 


± 0.05 


0.34 


± 0.02 


1.56 


2.69 ± 0.08 


0.59 ± 0.02 



Table 3. Details of the best-fitting exponential (Equation [H] dashed curves in Figure [3J and Plummer (Equation [5] dot-dashed curves 
in Figure models. Each scale radius is the geometric mean for the dwarf. The distance moduli and associated uncertainties derived in 
McConnachie et al. (2005) have been used to transform units. 



4 DISCUSSION 

Figures [2] and EJ and Tables Q |2j El and 0] present many 
results relating to the Cetus and M31 dSphs. In this sec- 
tion, we compare these results to those of previous studies 
(Section I4.1|l and then discuss each dSph individually (Sec- 
tions'^] and ^3J. We then compare the M31 and Galactic 



populations (Section 14.41 and examine correlations in the 
Local Group population as a whole (Section 14,50 . 

4.1 Comparison with previous work 

ICaldwell etafl l|l992t and ICaldwelll jl99flh derived struc- 
tural parameters for Andromeda I, II, III, V, VI and 
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Integrated Photometry 
m v My V So ri (') ri (kpc) 

2 2 



Andromeda I 


12.7 


± 


0.1 


-11.8 


± 


0.1 


24.7 


24.8 


2.8 


0.60 


Andromeda II 


11.7 


± 


0.2 


-12.6 


± 


0.2 


24.5 


24.7 


5.6 


1.06 


Andromeda III 


14.4 


± 


0.3 


-10.2 


± 


0.3 


24.8 


24.7 


1.6 


0.36 


Andromeda V 


15.3 


± 


0.2 


-9.6 


± 


0.3 


25.3 


25.6 


1.3 


0.30 


Andromeda VI 


13.2 


± 


0.2 


-11.5 


± 


0.2 


24.1 


23.9 


1.8 


0.42 


Andromeda VII 


11.8 


± 


0.3 


-13.3 


± 


0.3 


23.2 


23.6 


3.3 


0.74 


Cetus 


13.2 


± 


0.2 


-11.3 


± 


0.3 


25.0 


24.8 


2.7 


0.60 



Table 4. The integrated photometric properties of the dSphs. The integrated apparent magnitude (my) is uncorrected for reddening, 
while the absolute magnitudes (My) and central surface brightnesses have been de-reddened. V a is the directly measured central surface 
brightness while So is the central surface brightness obtained by normalising the radial profile to the integrated luminosity of the dSph. 
The values for n are obtained by integration of the appropriate King profile, with the exception of Andromeda II where a two-component 

2 

model has been used (see Section l4.3l . and Cetus, where the exponential profile has been used. 



VII based upon the integrated light of these galaxies. 
IWhiting et alJ il999l) has derived the structural parame- 
ters for Cetus. Table lists the values for ri, my and 
My (corrected for the extinction measurem ents and dis- 
tances used here) from [Caldwell et alJ (I1992T) (Andromeda 
I, II and im.lCaldwelll Jl999l) (Andromeda V, VI and VII), 
and lWhiting et all Jl99Sh (Cetus). The value of r i for Ce- 
tus has been e valua ted by integrating the King profile in 
IWhiting et ai] (Il999l) . 

We measure ri for Cetus to be twice that of 
IWhiting et al.l il999l) . These authors trace the Cetus sur- 
face brightness profile to ~ 4arcmins, while the mea- 
surement here extends to ~ 1 2 arcmins. Add i tional ly, the 
limiting surface brightness of IWhiting et al.l 1^999) is ~ 
29.5 mags arcsec -2 (see their Figure 4), while our profile 
extends ~ 1.5 mags arcsec -2 deeper. As a result, we find 
Cetus extends much further than originally measured. The 
brighter luminosity we derive is fully consistent with this in- 
crease in size. Identical arguments apply for Andromeda V 
and VII, where we measure their radial profiles out to much 
large r radii and over a larger range in surface brightness 
than ICaldwelll il999l) . and find them to be more extended 
and therefore brighter. 

A similar reason to the one above is not responsible for 
the di fference in results for Andromeda II. ICaldwell et alJ 
( 1992) measures the surface brightness profile for this galaxy 
over a similar radial and surface brightness range as we do, 
yet we derive a value for r i that is more than twice as large 
as that derived bv lCaldwell et all l|l992l ). We believe this is 
a result of the complexity of the profile of Andromeda II, 
which we discuss in Section 4.3. 



4.2 The isolated dwarf spheroidal in Cetus 

In the Local Group, there is a preference for dSph galaxies to 
be found as satellites of either M31 or the Galaxy, whereas 
dwarf irregular galaxies are preferentially found in more iso- 



lated environments. The origin of this position - morphology 
relation is not yet known, although it may well be related to 
the ram pressure and tidal stripping of dwarf galaxies ne ar 
large galaxies by hot, gaseous haloes iMaver et al 1 (l2001allhl 
2005)). Cetus and Tucana are the only dSph galaxies in the 
Local Group which are not clearly satellites of either M31 or 
the Galaxy. Analysis of their properties, and comparison to 
the dSphs which are satellites, may prove useful in determin- 
ing the origin of this position - morphology relation and in 
determining which evolutionary processes can be attributed 
to environmental effects. 

The King model tidal radius of Cetus is rt = 6.6 kpc, 
making it apparently the largest dSph in the Local Group by 
a significant amount. However, the stellar radial profile out 
to 12 arcmins does not show any evidence of turning over, 
and the derived value of rt is well beyond the final data 
point. Therefore, the concept of a tidal radius in this case 
is probably misleading as Cetus shows no evidence of tidal 
truncation. Cetus is currently ~ 750 kpc from the Galaxy, 
which places it ~ 680 kpc from M31. If Cetus has spent most 
of its evolution in isolation then we would not expect it to be 
truncated; its profile is fully consistent with an undisturbed 
system. From Equation 7. and assuming the masses of M31 
and the Galaxy to be ~ 10 12 Mq , then if Cetus has a typical 
mass for a dSph (~ 5 x 10 7 Mq), the lack of a tidal radius out 
to 6 or 7 kpc implies that after its last major star formation 
episode, Cetus has never been much closer that ~ 200 kpc 
to either M31 or the Galaxy. 

Kinematic information for Cetus is required to compli- 
ment this data. A radial velocity measurement may show 
whether the motion of Cetus suggests it ever having inter- 
acted with the Galaxy or M31. In addition, a stellar velocity 
dispersion measurement will constrain its mass. This infor- 
mation will help decide whether Cetus poses a challenge for 
formation and gas-loss models of dSphs which currently de- 
pend upon interactions with larger systems. 
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ri (') my My 



Andromeda I 


Z.O 


1 O 7^ 
1Z. ( 


_j_ 


U.U / 


11 O 

— 11. o 


Andromeda II 


2.3 


12.71 


± 


0.16 


-11.6 


Andromeda III 


1.3 


14.21 


± 


0.08 


-10.4 


Andromeda V 


0.6 


15.92 


± 


0.14 


-8.9 


Andromeda VI 


1.4 


13.30 


± 


0.12 


-11.4 


Andromeda VII 


1.3 


12.90 


± 


0.27 


-12.2 


Cetus 


1.4 


14.4 


± 


0.2 


-10.1 



Table 5. Previous estimates of r 1 , my and My for the Cetus and M31 dSphs, taken from Caldwell et al. (1992) (Andromeda I, II and 

2 

III), Caldwell (1999) (Andromeda V, VI and VII), and Whiting et al. (1999) (Cetus). The values for My have been recalculated for the 
distance and reddening estimates used here. 



4.3 The M31 dwarf spheroidals 

Andromeda I is one of the largest and brightest M31 dSphs 
and displays the distinctive S-shape indicative of ongoing 
tidal disruption (Figure [5^) , suggesting the presence of low 
surface brightness tidal tails. Although the average position 
angle of Andromeda I is ~ 22° , its isophotes twist from ~ 0° 
in the centre to ~ 40° i n the outer region s, again indicating 
strong tidal disruption dChoi et al.ll2002l) . From the surface 
brightness profile alone (Figure there is no sign that 
Andromeda I is being disrupted; there is no indication of 
a break in the profile, whose presence would n ormally be 
interp reted as due to tidal heating/disruption llChoi et alJ 
2002). Of course, this signature might be found at a lower 
surface brightness threshold. 

Andromeda I shows the strongest evidence of tidal dis- 
ruption of all the M31 dSphs from its contours alone. Deeper 
photometry and kinematics of this object will reveal the full 
extent of its tidal disruption. Kinematic data will be par- 
ticularly interesting, as it is possible that Andromeda I is 
surrounded by only a relatively shallow dark matter halo. 
Otherwise, it is not obvious that the outer stellar regions 
should be so distorted by the M31 tidal field. 

Andromeda II, shown in Figure is huge and cir- 
cular. It is currently ~ 185 kpc from M31 and has a large 
King tidal radius (r t = 4.2 kpc). If its orbit is relatively cir- 
cular, then this could account for its appearance since the 
M31 tidal field at this distance should be small. The average 
ellipticity of Andromeda II is ~ 0.22, although its central re- 
gions are significantly more circular than this; the innermost 
three isophotes of Figure^ are visibly more "core-like" than 
its outer regions. 

The radial profile of Andromeda II (Figure |3Jj) is quite 
unusual. The formal reduced \ 2 value for the King, Plummer 
and exponential models are relatively poor; both the King 
and exponential profiles underestimate the surface bright- 
ness in the central few arcmins by nearly a factor of two, 
while the Plummer profile also underestimates the central 
surface brightness slightly while failing to match at inter- 
mediate radii. The radius at which the King profile diverges 




radius (arcmin) 



Figure 4. The best fit King model to the surface brightness pro- 
file of Andromeda II, ignoring the inner 2 arcmins. The formal 
value of the fit is good, and the divergence of the King profile to 
the data at small radii suggests the presence of a secondary core 
component. 



from the data is close to the radius where the "core" devel- 
ops in the isopleth map. 

Figure [I] is a King profile fit to the Andromeda II pro- 
file ignoring the inner 2 arcmins. The formal quality of the 
fit is much better in this range (\ 2 = 1.13) and diverges 
at the same radius at which the "core" is observed in Fig- 
ure This reveals a factor of two excess of stars at small 
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radius, which is compelling evidence for a second compo- 
nent, in the form of a roughly constant density stellar core. 
Deeper global photometry down to the level of the hori- 
zontal branch and more kinematic information is required 
to confirm this hypothesis and place the structure of this 
galaxy in an evolutionary context. For example, it could be 
that this represents a dissolved star cluster, or the secondary 
component may reflect spatial variations in the star forma- 
tion history, such that there were two main episodes where 
one led to a more centrally concentrated stellar population 
than the other. I ndeed , from their HST study of this galaxy, 
iDa Costa et alJ <ll996t) are unable to model its abundance 
distribution without assuming a two component ("metal- 
rich" and "metal-poor") chemical enrichment model. The 
above scenarios, however, have only recently started to be 
conside red for the Galactic dSp hs, motivated by the discov- 
eries bv lKlevna et al.ll2004l and lTolstov et alJ ll2004) . 

Andromeda III is the most elongated dSph in this 
sample. Its outer isophote in Figure |2Ji suggests that the 
outer regions of this galaxy may be somewhat stretched 
and perturbed, although there is no evidence for this in the 
surface brightness profile (Figure |3J;) . Deeper photometry 
should reveal the true extent of this sparse population. An- 
dromeda III is one the closest dSphs to M31, at a distance 
of ~ 75kpc. 

Andromeda V (My ~ —9.6 mags) is a small and faint 
dSph. The main body of the dwarf is fairly round and com- 
pact, but its outermost isophotes appear diffuse and fuzzy 
(Figure |5Ji). It is possible that this is due to heating of the 
outer region by the tidal field of M31. The radial profile 
of Andromeda V appears to show a break at 3.5 arcmins, 
resulting in a failure of all three parametric fits to follow 
the profile. Simulations show that tidal effects can mani- 
fest themsel ves in the form o f a break in the surface bright- 
ness profile iChoi et al]l2002|) . and so Andromeda V may be 
showing evidence of tidal harassment. 

It is also worth noting that the absolute magnitude 
of Andromeda V which we derive is broadly consistent 
with a metallicity of [Fe/H] — —1.5 iArmandroff et all 
Il998l: IMcConnachie et al.l l2005h . Andromeda V was pre- 
viously measured to be fainter (Section 4.1), and its 
original metallicity me asurement of [Fe/H] ~ —1.5 by 
lArmandroff et all (I1998T) meant that this object was an out- 
lier on the lum inosity - metallicity relation first shown by 
ICaldwehl |l999). A later measurement marked this object a s 
more metal poor, at [Fe/H] ~ -2.2 foavidee et al.ll2002t) . 
more consistent with its low-luminosity measurement. How- 
ever, the most recent m easurement of its metallicity by 
lMcConriagIrie_ej i _aL| ( 2005) agrees with the original measure- 
ment of TATmandroffeT^dl l)l998f) . Given its absolute magni- 
tude measured here and its higher metallicity measurement, 
Androme da V appears to follow the luminosity - metallicity 
relation of lCaldwelJ l|l999T ) , while its earlier discrepancy can 
be attributed to an underestimated luminosity. 

Andromeda VI has a radial profile that is well de- 
scribed by a King profile. It shows no evidence of tidal dis- 
ruption in either the isopleth map (Figure^) or the surface 
brightness profile (Figure Et). Its ri is average for the M31 
dSphs, and its My and central surface brightness are like- 
wise relatively typical of the population. 

Andromeda VII is the brightest and most extended of 
the M31 dSphs. In a similar way to Andromeda II, Figure^ 



shows it to be nearly circular in projection, with an ellip- 
ticity of 0.13 and rt = 4.3 kpc. Andromeda VII is currently 
~ 220 kpc from M31 and a large tidal radius is not unex- 
pected if its orbit is relatively circular. Unlike Andromeda II, 
its surface brightness profile shows no deviation from a King 
model fit. 



4.4 Contrasting the M31 and Galactic dwarf 
spheroidals 

The M31 sub-gr oup consists of seven dSph galaxies includin g 
Andromeda IX JZucker et alJ 12004 IChapman et"aH 12005). 
compared to nine firmly identified dSphs around the Galaxy. 
We do not consider Sagittarius here, however, due to its 
extreme nature as a highly disrupted system near to the 
Galactic disk. 

IH95, and subsequent authors, found that King and ex- 
ponential profiles generally fit the surface brightness profiles 
of the Galactic dSphs very well. For the dSphs analysed here, 
we find King models are usually marginally better fits than 
exponential and Plummer profiles, which is not unexpected 
given that King profiles have an additional parameter. The 
latter two profiles fit the data equally as well as each other. 
Six out of eight of the galaxies looked at by IH95 show an ex- 
cess of stars at large radii in comparison to the best-fitting 
King model. Of the seven dSphs analysed here, only An- 
dromeda V shows clear evidence of this effect. Interpretation 
of these stars is a controversial topic in Local Group astron- 
omy. However, in the case of Andromeda V, where the stellar 
excess is a result of a break in the surface brightness pro- 
file, it is like ly that this is a m anifestation of tidal effects, 
as shown by Ichoi et al.1 J2002D . If this can also be shown 
to be the case for the relevant Galactic dSphs, then this 
might suggest that the outer regions of the Galactic dSphs 
are generally more disrupted than the M31 dSphs. However, 
only small samples are being compared and the evidence is 
indirect at best; Andromeda I, III and V all hint, to vary- 
ing degrees, that these galaxies may be tidally perturbed. 
It should also be noted that the photometry on which the 
various analyses of the Galactic dSphs are based extends 
below the main sequence turn-off, whereas the INT WFC 
photometry of M31 dSphs only samples bright RGB stars. 

Figure shows the absolute magnitude of each dSph 
plotted against ri (left panel) and r t (right panel) for most 
of the Local Group dSphs. Red squares represent the Galac- 
tic dSphs, blue triangles represent the M31 dSphs, and Cetus 
is shown as a magenta diamond . Da ta are taken from thi s 
study, IH95. IZucker et alJ J2004I) and lHarbeck et alJ J2005I) . 
The half-light radius of Andromeda IX is cal culated by inte- 
gration of the best fit King model deriv ed bvlHarbeck et af 
J2005J), using the distance calculated in IMcConnachie et al, 
(2005). For consistency, we have used the scale-radii corre- 
sponding to the single component fit for Andromeda II. 

Figure|^|demonstrates that the M31 and Galactic dSphs 
have a similar range of Mv\ the only noticeable difference 
in this respect is that there are five Galactic satellites with 
M v > —10 compared to two for M31. Although this is a pos- 
sible indication of relative incompleteness of the M31 dSph 
population, the small numbers involved mean that a K-S test 
shows that the two distributions of Mv are still consistent 
at the 62 % level. 
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Figure 5. Absolute magnitude versus ri (left panel) and rt (right panel) for the Galactic dSphs (red squares), M31 dSphs (blue triangles) 

2 

and the isolated dSph in Cetus (magenta diamond). For a given M v , the M31 dSphs have scale-radii that are generally at least twice as 
large as those for the Galactic dSphs. 



An obvious feature of Figure is that the values for 
ri and r t of the M31 dSphs are much larger than for their 
Galactic counterparts. For a given My, the scale-radii of 
the M31 dSphs are generally twice as large than for the 
Galactic dSphs. Specifically, the mean and median r t for 
the M31 dSphs are 2.0 and 2.3 times larger, respectively, 
than for the Galactic dSphs and the mean and median r i 

2 

are both 3.1 times larger. The half-light radius, ri, tracks 
r c , r e and b, and so the same disparity is also seen in these 
quantities. As these differences are observed across the range 
of My presented by the dSphs, it is unlikely to be an artifact 
of small number statistics. Instead, these findings point to 
notable differences between the formation and/or evolution 
of these dSph populations. 

A significant difference in the tidal radius, rt, between 
the two populations will, by default, lead to a difference in 
the values for r i . As demonstrated by the value for Cetus, 
rt depends strongly on environment. The value of rt de- 
rived from a King profile fit is not necessarily the true rt for 
a dSph, but is a useful paramete r isatio n to compare with 
simple analytic models. lOh et all il995l) give the following 
expressions for the value of r t for a dSph in the tidal field 
induced for a point mass; 



n 



(1 



(3 + e)* 



M 
Mi 



(7) 



stars in the dSph are taken into account, as stars on radial 
orbits are prefere ntially stripped to stars on circular orbits 
jRead et alJl2005l) . 

If we make the plausible assumption that the only fun- 
damental difference between the M31 and Galactic dSphs is 
that one group orbits M31 and the other group orbits the 
Galaxy, then the differences in rt (and the more robust mea- 
sure ri) must primarily arise via some combination of the 
relative distribution of semi-major axes, a modulated by the 
orbital eccentricities e, or the dynamical mass distribution 
embodied in the parameter Mh- However, it is difficult to 
break the degeneracy of these factors without more detailed 
modelling. For example, if Mh as a function of radius is dif- 
ferent between M31 and the Galaxy, then so too will the 
typical values of a. Alternatively, if the characteristic orbits 
of the two populations have significantly different values for 
e, then the fraction of the mass of the host galaxy contained 
within the positions of the dSphs could be very different and 
change as function of orbital phase, 9. Naively, however, it 
seems probable that the difference in the typical values of 
rt between the dSphs of M31 and the Galaxy is reflecting a 
difference in Mh (r, 9) between these two hosts, and requires 
detailed examinati on. In this context, i t is particularly inter- 
esting to note that iHuxor et alJ i2004) have recently found 
several extended luminous star clusters in the halo of M31 
with large core and tidal radii, which do not have any Galac- 
tic counterparts. 



and for a logarithmic halo, 



n 



(1 - ef 



+ 1 



MY 
Mh~) 



(8) 



Here, e is the eccentricity of the orbit of the dSph which has 
a semi-major axis a. M is the total mass of the dSph and Mh 
is the mass of the host galaxy contained within the current 
position of the dSph. The details of these equations are not 
important for the current discussion and we include them 
only to highlight the factors on which rt depends. Detailed 
treatments of rt also require that the orbits of the individual 



4.5 The Local Group population of dwarf 
spheroidals 

The Local Group dwarf galaxies display several correlations 
between their physical properties. These include correlations 
between luminosity - halo virial velocity, luminosity - metal- 
licity, spin-parameter - central surface brightnes s, and cen - 
tral surface brightness - luminosity. iDekel fc Wool ll200otl 
propose that all of these correlations relate to the role of su- 
pernovae feedback in dwarf galaxies , and is an extension of 
the ideas proposed bv lDekel fc Silkl jl986l) . The correlation 
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Figure 6. The relation between integrated luminosity and central 
surface brightness shown by the Local Group dSphs. For the M31 
dSphs, the central surface brightness is that which has been di- 
rectly measured from the integrated flux distribution. For a given 
My, the M31 dSphs have systematically fainter central surface 
brightnesses, presumably related to their large physical extent. 
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Figure 7. Left panel: central surface brightnesses of the Local 
Group dSph satellites as a function of distance from the host 
galaxy. For the M31 dSphs, the central surface brightness is that 
which has been directly measured from the integrated flux distri- 
bution. The central surface brightnesses of the Local Group dSph 
satellites seem to correlate with their current separation from 
their host, such that more distant dSphs are brighter. A Spear- 
man rank-order correlation test indicates that this correlation is 
significant at the 99 % level. 



between central surface brightness and luminosity is shown 
in Figure|S| and is significant at > 99.99 % level. While both 
the M31 and Galactic dSphs show the same trend, the two 
populations appear offset from one another such that the 
M31 dSphs have systematically lower central surface bright- 
nesses for a given Mv- This is presumably related to their 
more extended nature discussed in the previous section. 

Figure [7| shows the central surface brightness of the Lo- 
cal Group dSph satellites plotted against their distance from 



the host galaxies. There is an apparent correlation between 
central surface brightness and current separation, such that 
brighter dSphs are further from their host. A Spearman 
rank-order correlation test indicates that this correlation is 
significant at the 99 % level. Due to the correlation between 
central surface brightness and My- highlighted previously, 
there is a secondary correlation between Mv and separation 
which is significant at the 96 % level. 

This trend has already be en noted for the Galactic 
satellites; lBellazz in^e^l.1 1^961) fav oured a physica l expla- 
nation, whereas " iMateoT il99St) and Ivan den Berghl dl999l) 
suggested that this reveals that the Galactic satellite sys- 
tem is in complete at the faint end at large distances. Most 
recently, IWillman et al.l ll2004l) have suggested that some 
Galactic dSphs fainter than 24 mag arcsec -2 at a distance 
greater t han 100 kpc may not have been able to be detected. 
However. iMcConnachie fc Irwinl J2005) argue that Galactic 
satellites as faint as Sextans would have been found within 
200 kpc of the Galaxy in the 2/3 of the sky that were anal- 
ysed as part of the sur vey that led to the discovery of Sex- 
tans Jlrwin et al.lll990l) . The detection limits of this survey 
are discussed in detail in Irwin) (1199411 ; briefly, for satellites 
located between ~ 30 kpc and ~ 400 kpc from the Galaxy, 
the detectability of their resolved image signature is, to first 
order, oc S x A , where the first term is the central surface 
brightness and the second term is their scale-area. A oc d~ 2 

1 3 

and So oc d? , implying detectability scales as d~ 5 . This im- 
plies that: Fornax, Sculptor and Leo I would be detectable 
out to 400 kpc; Leo II would be detectable out to 350 kpc; 
Carina and Draco out to 300 kpc; Sextans and Ursa Mi- 
nor out to 250 kpc; satellites one magnitude fainter than 
Sextans or Ursa Minor would be detected out to 200 kpc. 
Given this, it is unlikely that the observed trend for the 
Galactic dSphs is a result of selection effects. Further, the 
M31 dSphs are here shown to display the same trend as the 
Galactic dS phs, although different select ion effects apply to 
this system l)McConnachie fc Irwinl2005t) . Clearly, this trend 
deserves further attention. 

It is important to consider the possibility that the cor- 
relation between central surface brightness and current sep- 
aration may be spurious, since the observed separation of a 
dSph from its host galaxy is a function of orbital phase. The 
current galactocentric distance of a satellite from its host 
may have no physical relevance in terms of the formation 
and evolution of that system. To try to quantify this effect, 
we have examined the Keplerian case of a satellite orbit- 
ing a point mass in order to determine for what fraction of 
an orbital period the instantaneous separation of a satellite 
from its host is a reasonable estimate of the time-averaged 
separation. This is equivalent to asking what fraction of a 
population of satellites on similar orbits will be at a distance 
from their host which is a reasonable representation of their 
mean distance, when viewed at a random moment in time. 

At a given moment, t, during an orbit of period P, the 
orbital phase is given by 6 = t/P and we can iteratively 
solve Kepler's equation for the eccentric anomaly, E, 

E = 2n9 + esinE , (9) 

for an orbit of eccentricity e. This allows the calculation of 
the instantaneous radius vector in units of the semi-major 
axis, 
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e cos E 



(10) 



The time-averaged separation, (r 1 ), of the satellite is then 
given by 



r') = / r'dO 
Jo 



(11) 



and the fraction of time, f e (77), for which r' is within a 
certain fraction of (r') is 



«2 



■qdO 



(12) 



77 = r'/ (r') and (81,02) corresponds to the range of for 
which r' is within the required fraction of (r'). The sub- 
script indicates that this is for a specific value of e, and the 
formalism is readily extended to a range of e since 



p(n) 



Jf e ( V )P(e)de 
f P (e) de 



(13) 



For simplicity, we set P (e) to be constant over some range 
of e. 

Table |S| illustrates the results for four differ ent satel- 
lite po pulations, each with a different spread in e. lOh et alJ 
(1995) estimated the range of e occupied by the Galactic 
satellites to be 0.3 < e < 0.7. Even for the case where orbits 
are radially biased, r' is still a reasonable indicator of (r') 
(ie. within 33%) in > 2/3 of the cases. These idealised cal- 
culations suggest that a correlation of r' with central surface 
brightness for the Local Group dSph population implies a 
similar correlation between (r') and central surface bright- 
ness. This then implies that an explanation for the trends is 
root ed in physics . 

lAbadi et alJ (l2005f) suggest that, in the current hierar- 
chical paradigm, the dwarf galaxies observed as satellites 
today were accreted relatively recently, and did not form in 
situ. If this is the case, then the above correlation is unlikely 
to be a result of the formation of the dwarfs but is more 
likely to be due to their subsequent evolution. It is difficult 
to significantly modify the characteristics of the central re- 
gions of a dwarf galaxy with an external tidal field, without 
severely disturbing/disrupting much of the system, since the 
inner regions of dSphs are more tightly bou nd and shielded 
than the outer regions fe.g. lOh et alJ |l995). The most sig- 
nificant interactions between the satellites and their hosts 
are gravitational. While ram pressure effects and supernova- 
driven winds undoubtedly contribute to the development of 
the baryonic density distribution, it seems likely that grav- 
itational effects play a major role in the evolution of the 
central surface brightness of satellites as they orbit and in- 
teract with their host. These effects are directly amenable 
to study by numerical simulations. 



shows a large variation in morphology, magnitude and ra- 
dial extent. Andromeda III and V show tentative evidence 
of tidal harassment, while the morphology of Andromeda I 
clearly identifies it as a strongly disrupted satellite of M31, 
suggesting that it may have extended tidal tails and perhaps 
only a relatively shallow dark matter potential. The isolated 
dSph in Cetus, on the other hand, does not show any evi- 
dence of tidal truncation and has a large radial extent. These 
results suggest a wide range of tidal effects is experienced by 
Local Group dSphs. Andromeda II has a clear excess of stars 
in its central regions and provides compelling evidence for a 
multiple- component system ; similar perhaps to Sculptor or 
Sextans dTolstov et al.ll2004 iKlevna et aljbooi) . Several of 
the dSphs are found to be more extended and brighter than 
previous studies, and in particular Andromeda V is shown to 
have a l uminosity consisten t with the metallicity measure- 
ments of lArmandroff et aT] il998t) and iMcConnachie et alJ 
(2005). Its previ ous value was consiste nt with a lower metal- 
licity derived bv lDavidee et alJ (I2002T) . 

In general terms, the M31 dSphs show qualitatively sim- 
ilar structural characteristics to the Galactic population. 
However, the M31 dSphs have much larger scale-radii for 
a given My than the Galactic dSphs, by a factor of 2 - 3. 
We suggest that these differences may result from differences 
in the tidal fields experienced by the two populations, such 
that the stronger Galactic tidal field truncates its satellites 
at smaller average radii in comparison to the M31 dSphs 
which evolve in a comparatively weaker tidal field. This im- 
plies the masses of M31 and the Galaxy differ and/or the 
characteristic orbits of the satellites are different, and re- 
quires detailed modellin g. In this context, it is particularly 
interesting to note that iHuxor et all (120041) have recently 
found several luminous extended star clusters in the halo of 
M31 with large scale radii, which do not have any Galactic 
counterparts. Finally, we highlight a correlation between the 
central surface brightnesses and the current separations of 
the dSphs from their hosts, which we argue is due to some 
unknown physical mechanism, which is most likely dynami- 
cal in origin. 



ACKNOWLEDGEMENTS 

This work is based on observations made with the Isaac 
Newton Telescope on the Island of La Palma by the Isaac 
Newton Group in the Spanish Observatorio del Roque de 
los Muchachos of the Institutode Astrofisica de Canarias. 
We would like to thank the anonymous referee for his/her 
comments, which led to an improvement in the clarity of 
this paper. Also, we thank Beth Willman for her feedback 
on an earlier draft. 



5 SUMMARY 

In this paper, we have derived isopleth maps, surface bright- 
ness profiles, intensity-weighted centres, position angles, el- 
lipticities, tidal radii, core radii, concentration parameters, 
exponential scale lengths, Plummer scale lengths, half-light 
radii, absolute magnitudes and central surface brightnesses 
for the Andromeda I, II, III, V, VI, VII and Cetus dSph 
galaxies. We have shown that the M31 dSph population 
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